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It has been shown numerically that systems of particles interacting with isotropic “stealthy” 
bounded long-ranged pair potentials (similar to Friedel oscillations) have classical ground states that 
are (counterintuitively) disordered, hyperuniform, and highly degenerate. Disordered hyperuniform 
systems have received attention recently because they are distinguishable exotic states of matter 
poised between a crystal and liquid that are endowed with novel thermodynamic and physical 
properties. The task of formulating an ensemble theory that yields analytical predictions for the 
structural characteristics and other properties of stealthy degenerate ground states in d-dimensional 
Euclidean space M'* is highly nontrivial because the dimensionality of the configuration space depends 
on the number density p and there is a multitude of ways of sampling the ground-state manifold, 
each with its own probability measure for finding a particular ground-state configuration. The 
purpose of this paper is to take some initial steps in this direction. Specifically, we derive general 
exact relations for thermodynamic properties (energy, pressure, and isothermal compressibility) 
that apply to any ground-state ensemble as a function of p in any d, and we show how disordered 
degenerate ground states arise as part of the ground-state manifold. We also derive exact integral 
conditions that both the pair correlation function g 2 {r) and structure factor S{k) must obey for 
any d. We then specialize our results to the canonical ensemble (in the zero-temperature limit) by 
exploiting an ansatz that stealthy states behave remarkably like “pseudo” equilibrium hard-sphere 
systems in Fourier space. Our theoretical predictions for 32 (r) and S{k) are in excellent agreement 
with computer simulations across the first three space dimensions. These results are used to obtain 
order metrics, local number variance and nearest-neighbor functions across dimensions. We also 
derive accurate analytical formulas for the structure factor and thermal expansion coefficient for 
the excited states at sufhciently small temperatures for any d. The development of this theory 
provides new insights regarding our fundamental understanding of the nature and formation of low- 
temperature states of amorphous matter. Our work also offers challenges to experimentalists to 
synthesize stealthy ground states at the molecular level. 

PACS numbers: 05.20.-y, 82.35.Jk,82.70.Dd 61.50.Ah 


I. INTRODUCTION 

The equilibrium structure and phase behavior of soft 
matter systems span from the relatively simple, as found 
in strongly repulsive colloidal particles, to the hi ghly 
complex, as seen in microemulsions and polymers [iHlOj . 
Soft matter has been fruitfully microscopically modeled 
as classical many-particle systems in which the particles 
(or metaparticles) interact with effective pair potentials. 
Bounded (soft) effective interactions have been particu¬ 
larly useful in modeling polymer systems, and they dis¬ 
play zero-temperature ground states with a rich variety 
of crystalline structures, depending on the coinposition of 
the constituents and interaction parameters 
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We have previously used a “collective-coordinate” ap¬ 
proach to generate numerically exotic classical ground 
states of many particles interacting with certain bounded 
isotropic long-ranged pair potentials in one-, two- and 
three-dimensional Euclidean space dimensions [TiHil as 
well as with anisotropic potentials [isj . It was shown 
that the constructed ground states across dimensions are 
the expected crystal structures in a low-density regime 
[ill [iu [13, but above some critical density, there 
is a phase transition to ground states that are, coun¬ 
terintuitively, disordered (statistically isotropic with no 
long-range order), hyperuniform, and highly degenerate 
[ 23 . These unusual amorphous states of matter have 
been shown to be endowed with novel thermodynamic 
and physical properties [l3 - [l6l . [23 - 1^ and belong to the 
more general class of disordered “hyperuniform” systems, 
which have been attracting attention recently, as detailed 
below. 
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The disordered ground states are highly degenerate 
with a configurational dimensionality that depends on 
the density, and there are an infinite number of distinct 
ways to sample this complex ground-state manifold, each 
with its own probability measure. For these reasons, it 
is theoretically very challenging to devise ensemble the¬ 
ories that are capable of predicting structural attributes 
and other properties of the ground-state configurations. 
A new type of statistical-mechanical theory must be in¬ 
vented to characterize these exotic states of matter. The 
purpose of this paper is to take some initial steps in this 
direction. However, to motivate the theoretical formal¬ 
ism, it is instructive to first briefly review the collective- 
coordinate numerical procedure that we have used to 
achieve disordered ground states. 

In the simplest setting, we previously examined pair¬ 
wise additive potentials v(r) that are bounded and inte- 
grable such that their Fourier transforms 'O(k) exist. If 
JV identical point particles reside in a fundamental region 
F of volume uf in d-dimensional Euclidean space at 
positions r-^ = ri,..., rjv under periodic boundary con¬ 
ditions, the total potential energy can be expressed in 
terms of u(k) as follows: 


$(r^) 


1 

2uf 


^u(k)|h(k)p 


iV^^(k) 

k 


( 1 ) 


where n(k) = exp(—fk • r^) is the complex collec¬ 

tive density variable, which can be viewed as a nonlinear 
transformation from the finite set of particle coordinates 
ri,..., rjY to the complex functions h(k) that depend on 
the infinite set of wave vectors k in reciprocal space ap¬ 
propriate to the fundamental cell F. The crucial idea is 
that if 'D(k) is defined to be bounded and positive with 
support in the radial interval 0 < |k| < AT and if the 
particles are arranged so that |h(k)p, a quantity propor¬ 
tional to the structure factor >S'(k), is driven to its min¬ 
imum value of zero for all wave vectors where 'O(k) has 
support (except k = 0), then it is clear from relation ([T]) 
that the system must be at its ground state or global en¬ 
ergy minimum. We have referred to these ground-state 
configurations as “stealthy” [I^ because the structure 
factor S'(k) (scattering pattern) is zero for 0 < |k| < K, 
meaning that they completely suppress single scattering 
of incident radiation for these wave vectors and, thus, 
are transparent at the corresponding wavelengths [2^. 
Various optimization techniques were employed to find 
the globally energy-minimizing configurations within an 
exceedingly small numerical tolerance [uHIl. Gener¬ 
ally, a numerically obtained ground-state configuration 
depends on the number of particles JV within the fun¬ 
damental cell, initial particle configuration, shape of the 
fundamental cell, and particular optimization technique 
employed. 

As the number of k vectors for which |h(k)| is con¬ 
strained to be zero increases, i.e., as K increases, the di¬ 
mensionality of the ground-state configuration manifold 
Fc decreases. Because |n(k)| is inversion symmetric, the 


number of wave vectors contained in a sphere of radius 
K centered at the origin must be an odd integer, say 
2M{K) F I, and thus M{K) is the number of indepen¬ 
dently constrained wave vectors [2^ . The parameter 


M{K) 
d{N -ly 


( 2 ) 


which is inversely proportional to density, gives a mea¬ 
sure of the relative fraction of constrained degrees of free¬ 
dom compared to the total number of degrees of freedom 
d{N — I) (subtracting out the system translational de¬ 
grees of freedom). We show in Sec. lIVBI that the dimen¬ 
sionality of the configuration space per particle is given 
by d{l — 2x) in the thermodynamic limit. 
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FIG. 1. Schematic illustrating the inverse relationship be¬ 
tween the direct-space number density p and relative fraction 
of constrained degrees of freedom y for a fixed reciprocal- 
space exclusion-sphere radius K (where dark blue k points 
signify zero intensity with green, yellow, and red points in¬ 
dicating increasingly larger intensities) for a stealthy ground 
state. A compression of a disordered ground-state configura¬ 
tion with a hxed number of particles N in direct space leads 
to a dilation of the lattice spacing in reciprocal space. This 
means that during the compression process, the k points for 
which |h(k)| is zero associated with the initial uncompressed 
system move out of the exclusion zone; i.e., the value of M{K) 
[cf. Eq. ([2])] decreases. Since there are fewer constrained de¬ 
grees of freedom (dimensionality of the ground-state config¬ 
uration manifold increases), the disordered direct-space con¬ 
figuration becomes less spatially correlated. For a fixed N in 
the limit p —>■ oo (i.e. system volume vf —> 0), every k point 
(except the origin) is expelled from the exclusion zone, and 
the system tends to an ideal-gas configuration [s^, even if it 
is not an ideal gas thermodynamically, as shown in Sec. IIV Bl 


It is straightforward to see why, for sufficiently small 
X, ground states exist that are highly degenerate and 
typically disordered for sufficiently large N-, see Fig. [T] 
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Clearly, when the system is free of any constraints, i.e., if 
X = 0, it is a noninteracting classical ideal gas. While it 
is unusual to think of classical ideal-gas configurations as 
ground states, at T = 0, they indeed are global energy¬ 
minimizing states that are highly degenerate and typi¬ 
cally disordered for large enough N. While the ground- 
state manifold contains periodic configurations (e.g., Bra- 
vais lattices and lattices with a basis), these are sets of 
zero measure in the thermodynamic limit. Clearly, if x is 
made positive but very small, the ground states remain 
disordered and highly degenerate, even if the dimension¬ 
ality of the configuration space Dc has now been sud¬ 
denly reduced due to the imposed constrained degrees 
of freedom, the number of which is determined by the 
radius K. From relation we see that if K is fixed, 
configurations with ideal-gas-like pair correlation func¬ 
tions correspond to the limit x 0 or, equivalently, to 
the limit p —^ cx) [s^. The latter situation runs counter 
to traditional understanding that ideal-gas configurations 
correspond to the opposite zero-density limit of classical 
systems of particles. The reason for this inversion of lim¬ 
its is due to the fact that a compression of the system in 
direct space leads to a dilation of the lattice spacing in 
reciprocal space, as schematically shown in Fig. [TJ While 
it is not surprising that the configuration space is fully 
connected for sufficiently small x, quantifying its topol¬ 
ogy as a function of x for all allowable x is an outstanding 
problem, which is discussed further in the Conclusions. 


FIG. 2. A disordered nonhyperuniform configuration (left 
panel) and a disordered hyperuniform configuration (right 
panel). We arrive at the configuration on the right by very 
small collective displacements of the particles on the left via 
the methods described in Ref. [l^. (Each particle on av¬ 
erage moves a root-mean-square distance that is about an 
order of magnitude smaller than the mean-nearest-neighbor 
distance as measured by the configuration proximity metric 
[3ll|.i These two examples show that it can be very difficult to 
detect hyperuniformity by eye, and yet their large-scale den¬ 
sity fluctuations are dramatically different (“hidden order”). 


It is noteworthy that stealthy point patterns (disor¬ 
dered or not) constitute a special class of so-called hy¬ 
peruniform states of matter. Hyperuniform systems are 
characterized by vanishing (normalized) density fluctu¬ 
ations at large length scales; i.e., the structure factor 
S'(k) tends to zero in the limit |k| —0 [s^ (see Sec. 
Ill Al for details). The hyperuniformity concept provides 
a means of categorizing crystals, quasicrystals, and spe¬ 
cial disordered systems according to the degree to which 


large-scale density fluctuations are suppressed [H, [13 ■ 
Disordered hyperuniform patterns, of which disordered 
stealthy systems are special cases, behave more like crys¬ 
tals in the manner they suppress large-scale density fluc¬ 
tuations, and yet they also resemble typical statistically 
isotropic liquids and glasses with no Bragg peaks. In 
this sense, they have a “hidden order” on large length 
scales that is not apparent at small length scales, even 
if short-range order is present; see Fig. for a vivid 
illustration. During the last decade, a variety of disor¬ 
dered hyperuniform states have been identified that exist 
as both equilibrium and nonequilibrium phases, includ¬ 
ing maximally random jammed particle packings 
jammed athermal gran ular media [s^ - jammed thermal 
colloidal packings |39l l40l| . cold atoms [4l|, transitions 
in nonequilibrium systems 1^ , surface-enhanced Ra¬ 
man spectroscopy |44| . terahertz quantum cascade laser 
[ 43 , wave dynamics in disordered potentials based on su¬ 
persymmetry [13 , avian photoreceptor patterns [13 , and 
certain Coulombic systems [13. Moreover, disordered 
hyperuniform materials possess novel physical properties 
potentially important for applications in photonics psl - 
[ 23 , m, [13 and electronics [5l[ - [53 . 

The well-known compressibility relation from statisti¬ 
cal mechanics [13 provides some insights about the re¬ 
lationship between temperature T and hyperuniformity 
for equilibrium systems at number density p: 

S{k = 0) = pksTKT- (3) 

We see that any ground state (T = 0) in which the 
isothermal compressibility kt is bounded and positive 
must be hyperuniform because the structure factor S{k = 
0) must be zero. This includes crystals as well as exotic 
disordered ground states such as stealthy ones. However, 
in order to have a hyperuniform system at positive T, 
the isothermal compressibilit y m ust be zero; i.e., the sys- 
tem must be incompressible [iTj l (see Refs. [13 and [13 
for some examples). Subsequently, we will use relation 
© to draw some conclusions about the excited states 
associated with stealthy ground states. 

Our general objectiye is the formulation of a predic- 
tiye ensemble theory for the thermodynamic and struc¬ 
tural properties of stealthy degenerate disordered ground 
states in arbitrary space dimension d that c omp lements 
preyious numerical work on this topic [lll - ll6l| - After 
proyiding basic definitions and describing a family of 
isotropic stealthy potentials (Secs. HIl and iHll) . we de- 
riye general exact relations for thermodynamic proper¬ 
ties (energy, pressure, and isothermal compressibility) 
that apply to any well-defined ground-state ensemble as a 
function of the number density or, equivalently, x in any 
space dimension d (Sec. IIVI) . We subsequently derive 
some exact integral conditions that both the pair cor¬ 
relation function g 2 (j') and structure factor S{k) must 
obey (Sec. |V|. The existence of periodic stealthy ground 
states enables us to show how disordered degenerate 
ground states arise as part of the ground-state mani¬ 
fold for sufficiently small x (Sec. I VII) . Subsequently, we 
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derive analytical formulas for the pair statistics [g 2 if’) 
and S{k)] for sufficiently small x in the canonical ensem¬ 
ble in the limit that the temperature T tends to zero 
(Sec. IVIII) by exploiting an ansatz that stealthy states 
behave like “pseudo” equilibrium hard-sphere systems in 
Fourier space. Our theoretical predictions for g 2 {i’) and 
S{k) are in excellent agreement with computer simula¬ 
tions across the first three space dimensions. These re¬ 
sults are then used to predict, with high accuracy, other 
structural characteristics of stealthy ground states across 
dimensions, such as order metrics, local number variance, 
and nearest-neighbor functions (Secs. IVIIIl and llXI) . Sub¬ 
sequently, we derive analytical formulas for the struc¬ 
ture factor and thermal expansion coefficient for the as¬ 
sociated excited states for sufficiently small temperatures 
(Sec. 12). Finally, we provide concluding remarks in Sec. 

im 


II. DEFINITIONS AND PRELIMINARIES 

Roughly speaking, a point process in d-dimensional 
Euclidean space is a distribution of an infinite num¬ 
ber of points in with the configuration ri, r 2 ,... at a 
well-defined number density (number of points per unit 
volume). For a statistically homogeneous point process 
in at number density p [d^, the quantity p"g„(r") 
is the probability density associated with simultaneously 
finding n points at locations r" = ri, r 2 ,..., r„ in 
[s^l . With this convention, each n-particle correlation 
function gn approaches unity when all of the points be¬ 
come widely separated from one another. Statistical ho¬ 
mogeneity implies that gn is translationally invariant and 
hence only depends on the relative displacements of the 
positions with respect to any chosen system origin, e. 5 ., 
gn = g„(ri 2 , ri 3 ,..., ri„), where Vij = rj - r^. 

The pair correlation function 52 (r) is a particularly im¬ 
portant quantity. If the point process is also rotationally 
invariant (statistically isotropic), then g2 depends on the 
radial distance r = |r| only, i.e., 32 ( 1 *) = ff 2 (?')- Thus, it 
follows that the expected number of points Z{R) found 
in a sphere of radius R around a randomly chosen point 
of the point process, called the cumulative coordination 
function, is given by 

nR 

Z{R) = psi{l) / x‘^ ^g2{x)dx, (4) 

Jo 

where si(r) = fT{dj2) is the surface area of a 

d-dimensional sphere of radius r. The total correlation 
function h{r) is trivially related to 52 (r) as follows: 

d(r) =g 2 (r) - 1. (5) 

When there are no long-range correlations in the system, 
h{r) —> 0 or, equivalently, g 2 (r) ^ 1 as |r| —>• 00 . The 
structure factor <S'(k), which plays a prominent role in 
this paper, is related to the Fourier transform of h{r), 




FIG. 3. (Bravais) lattice with one particle per fundamental 
cell (left panel) and a periodic crystal with multiple particles 
per fundamental cell (right panel). 

denoted by d(k), via the expression 

.5(k) = l + pd(k). (6) 

A lattice A in R'^ is a subgroup consisting of integer lin¬ 
ear combinations of vectors that constitute a basis for R'^, 
and thus, it represents a special subset of point processes. 
Here, the space can be geometrically divided into iden¬ 
tical regions F called fundamental cells, each of which 
contains just one point specified by the lattice vector 

p = niai -I- n2a2 -I - \-nd-iad-i + Udau, (7) 

where a^ are the basis vectors for a fundamental cell and 
Hi spans all the integers for * = 1, 2, • • • , d. We denote by 
vf the volume of F. A lattice is called a Bravais lattice in 
the physical sciences. Unless otherwise stated, we will use 
the term lattice. Every lattice has a dual (or reciprocal) 
lattice A* in which the lattice sites are specified by the 
dual (reciprocal) lattice vector q • p = 27rm for all p , 
where m = 0, ±1, ±2, ±3 ■ • •. The dual fundamental cell 
F* has volume vp* = {‘2 ttY/vf- This implies that the 
number density pA of A is related to the number density 
PA» of the dual lattice A* via the expression 

PAPA* = l/(27r)‘^. (8) 

Some common d-dimensional lattices are mathematically 
defined in Appendix [2 

A periodic point process (crystal) is a more general 
notion than a lattice because it is obtained by placing a 
fixed configuration of N points (where A" > 1) within a 
fundamental cell E of a lattice A, which is then periodi¬ 
cally replicated. Thus, the point process is still periodic 
under translations by A, but the N points can occur any¬ 
where in F] see Fig. [3] 

A. Hyperuniform point processes 

Consider uniformly sampling the number of points that 
are contained within a spherical window of radius i? of a 
point process in R'^. A hyperuniform point process has 
the property that the local number variance cr^ (i?) grows 
more slowly than R'^ [s^. Because cr‘^{R) is exactly re¬ 
lated to a d-dimensional volume integral of the structure 
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factor S'(k) (see Sec. IVIIII) . this implies that hyperuni¬ 
form states of matter possess infinite-wavelength density 
fluctuations (appropriately normalized) that vanish; i.e., 
^(k) obeys the condition 

lim 5'(k) ^ 0, (9) 

IkKo 


v{k) are also radial functions in where fc = |k| is a 
wave number. The d-dimensional Fourier transform of 
any integrable radial function /(r) in is given by [H] 


fik) 


(27r)' 





J(d/ 2 )-i{kr) 

(fc^)W2)-l 


dr, 


( 11 ) 


which means they are poised at an “inverted” criti¬ 
cal point with associated scaling exponents [32|. For a 
Poisson (spatially uncorrelated) point process and many 
disordered point patterns, including typical liquids and 
structural glasses, the number variance grows like the 
volume of the window, i.e., cr^(i?) ^ R'^, implying that 
5'(k) is positive at k = 0 . All perfect crystals and qua¬ 
sicrystals are hyperuniform such that cr‘^{R) 
in other words, the variance grows like window surface 
area. By contrast, it is much more unusual to hnd dis¬ 
ordered systems that are also hyperuniform. In recent 
years, evidence has been emerging that disordered hype¬ 
runiform many-particle systems can be regarded as new 
distinguishable states of disordered matter (see examples 
given in the Introduction). Whenever the structure fac¬ 
tor goes to zero with the power-law form S'(k) ^ |k|“, 
the number variance has the following large-i? asymp¬ 
totic scaling that depends on the value of the exponent 


r R‘^-HnR 

a = 1 


a^{R) - < 

a < 1 

(R^+oo). (10) 


a > 1 



Since disordered as well as ordered stealthy states can 
be viewed as systems in which a tends to infinity, we 
see from Eq. OT that they have the asymptotic scal¬ 
ing <j‘^{R) ~ R‘^~^. We give theoretical predictions for 
the variance of disordered stealthy ground states in Sec. 

Enn 


and the inverse transform of f{k) is given by 


fir) 


1 

(27r)^ 



k'^-^hk) 


J{d/ 2 )-i{kr) 

(fc^)W2)-l 


dk, 


( 12 ) 


where Jy{x) is the Bessel function of order v. 

Consider the class of stealthy radial potential functions 
v{k) in R'^ that are bounded and positive with compact 
support in the radial interval 0 < k < K, i.e., 


v{k) = V{k)&{K -k), (13) 

where, for simplicity, V(k) is infinitely differentiable in 
the open interval [0, K) and 


0(x) 


0, X < 0 

1, X > 0 


(14) 


is the Heaviside step function. The corresponding direct- 
space radial pair potential v(r) is necessarily a delocal¬ 
ized, long-ranged function that is integrable in R*^. More¬ 
over, without any loss of generality, it will be assumed 
that V{k) < vq. 

For concreteness and purposes of illustration, we will 
examine properties of two specific families of potentials 
that fall within the aforementioned wide class of stealthy 
interactions: “power-law” and “overlap” potentials. 


A. Power-law potentials 


III. FAMILIES OF STEALTHY PAIR 
POTENTIALS 

As we see in the next section, the specific form of a 
stealthy potential does not affect the ground-state en¬ 
ergy manifold, but it can affect other thermodynamic 
properties, such as the pressure. This has consequences 
in simulations of such properties, especially with respect 
to convergence issues. Hence, it is instructive to remark 
on some mathematical aspects of the long-range nature 
of the direct-space stealthy potentials, which are very 
similar to the weakly decaying Friedel oscillations of the 
electron density in a variety of systems, including molten 
metals as well as graphene [HB, [E^- As we will see, in 
some cases, stealthy potentials may mimic effective in¬ 
teractions that arise in certain polymer systems 0- 
Here, we will limit ourselves to pair potentials v{r) that 
are radial functions in R"^, where x = |r| (i.e., isotropic 
pair interactions), and therefore, their Fourier transforms 


The power-law potentials are defined in Fourier space 
as follows: 

d{k)=vo{l-k/K)"^&{K-k), (15) 

where the exponent m can be any whole number. The 
corresponding direct-space potential v{r) will depend on 
d for any given m and is exactly given by 

vjr) _ K‘^T{m + l)F((d-|-1)/2) ■ 1 ^ 2 ( 01 ; 6 i, 62 ; a^) 

vq r(l-|-m-I-d)7r(‘^+i)/2 ’ 

(16) 

where ai = (d -I- l)/2, bi = 1 -I- (m -I- d)/2, 62 = 
(1 -I- m -I- d)/2, X = —{Kr)‘^/A and 1 ^ 2 ( 01 ; 61 , 62 ; x) is 
a special case of the generalized hypergeometric func¬ 
tion ^^^( 01 ,..., Op; 61 ,..., 5,; x) [ 5 ^. Because the potential 
(fTH)) is derived from the Fourier power-law potential (jl5|) , 
we will refer to Eq. (HU) as the direct-space power-law 
potential. In the instance when m = 0 in Eq. (Ha (i.e.. 
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simple step function), this expression for v(r) simplifies 
as follows: 


v(r) 

Vo 



d/2 

Jd/ 2 {Kr), 


(17) 


where s(r; m, d) is a bounded function (a sinusoidal func¬ 
tion or constant of order one). For any fixed d and 
1 < TO < d, the long-range oscillations of v{r) are 
controlled by an envelope that decays like 1/r^, where 
{d+l)/2 <P <d+l. 


for which the large-r asymptotic behavior is given by 


v(r) f cos(iFr — (d-I-l)7r/4) 

--{ii) - 

(18) 

For any fixed d and m > d, the direct-space power-law 
potential has the asymptotic form 


v{r) s(r;m,d) 

Vq j.d+1 


(r oo), 


(19) 


In Fig. |4l we plot the Fourier power-law potential for 
selected values of to (which applies in any dimension) 
and the corresponding direct-space potentials for d = 3. 
In all cases, we set vo = K = \. It is to be noted that 
the amplitudes of the oscillations in v{r) decrease as to 
increases for a fixed dimension. 



FIG. 4. Left Panel: Fourier power-law potential v{k) for the special cases m = 0, 2, and 4 that apply for any d. Right panel: 
Corresponding direct-space power-law potentials v(r) in the instance d = 3. Here, we set vo = K = 1. 


B. Overlap potentials 

Let a(r; R) represent the intersection volume of two 
identical d-dimensional spheres of radius R (scaled by 
the volume of sphere) whose centers are separated by 
a distance r. This quantity is known analytically in an^y 
space dimension, and has a variety of representations |59l| , 
including the following: 


i{r;R) = c{d) f 
Jo 


cos ^{r/(2R)) 


sin'^(6»)d6», (20) 


where c(d) is the d-dimensional constant given by 

2r(i + d/2) 

^ ’ ^i/2r((d + i)/2)' 

For d = I, 2, 3 and 4, we respectively have 


a{r-,R) = e{2R-r) 


I - 


2Ri 


( 21 ) 


( 22 ) 


r(r;i?) = 0(2i?-r) 



-1 f - — (i - 
\2RJ 2R \ 4i?2 J 

(23) 


a{r;R) = Q{2R-r) 


a{r;R) = Q{2R-r) 


1 - 


3 r 

a'r 


cos 


-if M_ 

\2R) 


5r 


1 / r 
16 \R 


1 / r 
12 \R 


(24) 




Consider the class of “overlap” potentials, which for 
any dimension is given by 


v{k) = Vo a{r = k,R = iF/2). 


(26) 


Note that for d = 1, the overlap potential is identical to 
the power-law potential when d = 1 and to = 1. The 
thermodynamics of the ground-state manifold of this po¬ 
tential in the special case d = 2 was numerically inves¬ 
tigated in Refs, [l^ and [l^. It follows from Eq. (|26ll 
that the corresponding direct-space overlap potential is 
given by 


7 (r) _ r(l + d/ 2 )^d/ 2 (^i'/ 2 ) 


Vo 


.d/2 


(27) 


1/2N 


which is clearly non-negative for all r. Its large-r asymp¬ 
totic behavior is given by 


'kl .. 4r(l + d/2) 

ypd ! 2 -\- \ 


Vo 


cos^ l^iFr/2 — 


(d -I- 1)7 


(r 


00 


(28) 
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revealing that the long-ranged decay of the direct-space 
overlap potential has an envelope controlled by the in¬ 
verse power law 

Figure [5] depicts the overlap potential v{k) for the 
first three space dimensions and the corresponding direct- 


space overlap potentials n(r), the latter of which vividly 
shows the increasing decay rate of v(r) with increasing di¬ 
mension. The direct-space overlap potential v(r) is sim¬ 
ilar in functional form to effective positive pair interac¬ 
tions that arise in multilayered ionic microgels 0- 




FIG. 5. Left panel: Fourier overlap potential v{k) for the first three space dimensions. Right panel: Corresponding direct-space 
overlap potentials v{r), which oscillate but are always non-negative. Here, we set vo = K = 1. 


IV. ENSEMBLE THEORY FOR STEALTHY 
DISORDERED GROUND STATES: EXACT 
RESULTS FOR THERMODYNAMIC 
PROPERTIES 

Our general objective is the formulation of an ensem¬ 
ble theory for the thermodynamic and structural proper¬ 
ties of stealthy degenerate disordered ground states that 
we previously investigated numerically [TTI - [l^ . In this 
section, we derive general exact relations for thermody¬ 
namic properties that apply to any well-defined ensemble 
as generated by a particular way to sample the stealthy 
disordered ground-state manifold as a function of num¬ 
ber density p. In the subsequent section, we derive some 
exact results for the pair statistics for general ensembles. 


A. Preliminaries 

To begin, consider a configuration of N identical par¬ 
ticles with positions = ri,... ,rAr in a large region 
of volume V in d-dimensional Euclidean space For 
particles interacting via a pair potential f(r), the total 
potential energy ^(r^) is given by 

$(r^) = ^u(r,,), (29) 

i<j 

where — r^. Of particular interest are classical 

ground states, i.e., those configurations that minimize 
the energy per particle $(r^)/iV. 

The ensemble average of the energy (12^ per particle 
u in the thermodynamic limit can be written in terms of 


the pair correlation function 52 (r): 

u = = R / w(r) 52 (r)dr 

= ^[ v{r)dr + ^[ v{r)h{r)dr,{30) 

where angular brackets denote an ensemble average and 
p is the number density in the thermodynamic limit. 

Because the collective-coordinate approach relies on 
the Fourier representation of the energy, we recast Eq. 
(l30)l in terms of the structure factor S'(k) by applying 
Parseval’s theorem to the second line of Eq. (l30l) : 

u=iv(k = 0) + ^J^Hk)h(k)Jk 

= |s(k = 0) - i.(r = 0) + ^ »(k)S(k)«) 

where ^(k) and h(k) are the Fourier transforms of v{r) 
and h{r), respectively, both of which are assumed to ex¬ 
ist, k is a wave vector, and >S'(k) is the structure factor 
defined in relation ([ 6 ]). Note that the structure factor is 
a non-negative, inversion-symmetric function, i.e., 

^(k) > 0 for all k, S'(k) = S'(-k). (32) 


B. Ground-state energy and dimensionality of its 
configuration space 

Consider a radial (isotropic) stealthy potential func¬ 
tion v{k) with support in 0 < fe < A of the class specified 
by Eq. (IT^ . In light of Eq. (ISTI) , it is clear that when¬ 
ever particle configurations in R'^ exist such that S{k) 
is constrained to achieve its minimum value of zero for 
0 < k < K, the system must be at its ground state or 




















global energy minimum. This follows because the inte¬ 
grand v{k)S{k) in the nontrivial term on the right side 
of Eq. (1^ is identically zero because of the conflicting 
demands of the step functions. When such configurations 
exist, the average ground-state energy per particle in any 
well-defined ensemble is given exactly by 

u=^vo- ^v{r = 0), (33) 

which is a structure-independent constant, depending on 
the density p, vq = v{k = 0), v{r = 0) and d, as explic¬ 
itly shown below. Importantly, because u is a constant 
independent of the structure, its value has no effect on 
the ground-state manifold, which is generally degenerate, 
; hence, this manifold is invariant to the specific choice 
of the stealthy function v(k) at fixed p and d. 

We would like to express the energy (l33)) in terms of the 
parameter %, defined by relation which measures the 
relative number of independently constrained degrees of 
freedom for a finite system under periodic boundary con¬ 
ditions. Note that, in the thermodynamic limit, M{K) in 
Eq. (m is simply half of the volume of a sphere of radius 
K [due to the inversion symmetry of S (k)] multiplied by 
the density pA» of the dual lattice [cf. Eq. dH)], i.e.. 


M{K) = PA. 


vi{K) 

2 


vi(K) 

2(27r)^p’ 


(34) 


where we have used the fact that p —>■ NpA in this dis¬ 
tinguished limit. Hence, from Eq. we obtain the 

following expression for y in the thermodynamic limit: 


where 


vi{K) 

2d(27r)'^’ 


vi{R) 


T^d/2j^d 

r(i-f d/2) 


(35) 


(36) 


is the volume of a d-dimensional sphere (hypersphere) of 
radius R. We see that for fixed K and d, which fixes the 
potential, x is inversely proportional to p, which is the 
situation that we usually consider in this paper (^ . 

Hence, as x tends to zero, p tends to infinity, which 
configurationally corresponds counterintuitively to the 
uncorrelated ideal-gas limit (Poisson distribution), as dis¬ 
cussed in the Introduction. As x increases from zero, the 
density p decreases and the dimensionality of the ground- 
state configuration manifold Dq decreases. The config¬ 
urational dimensionality per particle in the thermody¬ 
namic limit, dc, can easily be obtained from the relation 
Dc = dN — 2M{K) for a finite system [^; specifically. 


dc = d(l - 2x), (37) 


where dc = limDp_,.oo,Ar^.oo Dc/N. 

Equations (1551) and ((551) yield the average ground-state 
energy per particle to be 


u = Vo 


'P 

.2 


idpX , 


iP*miu <P<oo), (38) 


where p^^^ is the minimal density associated with the 
dual of the densest Bravais lattice in direct space (as 
elaborated in Sec. EH), and 

^ lu-i vik)dk ^ (27r)^^(r = 0) 

voVi{K) voVi{K) 

is a constant whose value depends on the specific form of 
the stealthy-potential class v{k) defined by Eq. (1T51) and 
hence must lie in the interval (0, 1 ], where 7 = 1 corre¬ 
sponds to the step-function choice v{k) = voQ{K — k). 
While the system in the limit x ^ 0 (p 00) corre¬ 
sponds configurationally to an ideal gas in so far as the 
pair correlation function is concerned, as we will explain 
in detail in Sec. lVIlAl thermodynamically, it is nonideal; 
see Eq. (1551) for u and Eq. (15T|) for the pressure. 


C. Energy route to pressure and isothermal 
compressibility 


The pressure in the thermodynamic limit at T = 0 can 
be obtained from the energy per particle via the relation 


P = P 


2 



(40) 


Therefore, for stealthy potentials, we see from Eq. (|551) 
that the ground-state pressure, for all possible values of 
p or X, is given by the following simple expression: 



(Pmm <P<00). 


(41) 


Hence, the isothermal compressibility 
of such a ground state is 


n-i 



T 



(42) 


We see that as p tends to infinity, the compressibility 
tends to zero. 

Two important remarks are in order. First, estimates 
of the pressure obtained from simulations that we previ¬ 
ously performed for d = 2 [mil , as well as those carried 
out in the present study across the first three space di¬ 
mensions (Appendix B), are in very good agreement with 
the exact result dm across a wide range of densities, thus 
validating the accuracy of the simulations. Second, the 
fact that the pressure (1411) is a continuous function of den¬ 
sity implies that any phase transition that may take place 
could be a continuous one, the implications of which are 
discussed in the Conclusions. 


D. Virial route to pressure and isothermal 
compressibility 

An alternative route to the pressure for a radial pair 
potential function v(r) is through the “virial” equation, 
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which at T = 0 in the thermodynamic limit, is given by 



dv 

dr 


g 2 {r)dr. 


(43) 


Although the pressure obtained via the virial route is 
generally expected to be equivalent to that obtained from 
the energy route (as described in the previous section), 
we will show that, for a certain class of stealthy poten¬ 
tials, the pressure obtained from Eq. (1431) is either ill 
defined or divergent. This has practical implications for 
what types of stealthy potentials can be used in constant- 
pressure simulations. 

It is convenient to rewrite the virial relation (1431) in the 
following form: 


P = - 


2d 


r—dv + 
dr 


r— h(r)dr 
dr 



F{k = 0) -f 


1 

[2'k)’^ 


F{k) h{k)dk 


(44) 


where F{k) is the Fourier transform of F{r) = rdv/dr, 
when it exists, and we have used Parseval’s theorem and 
definition (O for the total correlation function h{r). 

To continue with this analysis, we make use of the 
following lemma. 

Consider a bounded radial function z{k) with compact 
support on the radial interval [0, K\ in that is in¬ 
finitely differentiable in the open interval [0,K). There¬ 
fore, its Fourier transform z{r) exists. 

Lemma 1 .—^The Fourier transform of the radial func¬ 
tion w{r) = rdz/dr in is given by 


w{k) 


—d ■ z{k) 


dz 

dk 


(45) 


Proof.—Differentiation of z{k) [defined via Eq. dill) ] with 
respect to k leads to the following identity: 



( 27 r) 



J(d/ 2 ){kr 


dr. 


(46) 


The Fourier transform of w{r) is given by 


w{k) 


( 27 r) 


f°° gdz J(d/ 2 )-i{kr) 

Jo ^ dr 


(47) 


Integrating relation (gll) by parts and using Eq. (|46l) 
proves the lemma. 

Corollary .—It immediately follows from Lemma 1 that 
w{k) has the same support as z{k) and 


w{k = 0) = —d ■ z{k = 0), (48) 

meaning that the volume integral of rdz/dr over all space 
is proportional to the corresponding volume integral of 
z(r). 


Note that by the Corollary of Lemma 1, F(0) = —d • 
{i(0), and hence we can rewrite the virial relation (IT41) as 


= PJ- 

2d 

Z2 


dvik = 0) 


dvik = 0) — 


P 

= y.o. 



p{‘2.ttY Jr' 


F{k) h{k)dk 
F{k)dk 


(49) 


The second term in the second line of Eq. (HHl) follows 
because ph{k) = —Q{K — k) inside the exclusion sphere 
of radius K [see also Eq. (l5^ below] and has support in 
this exclusion zone by the Corollary of Lemma I. But this 
second term must vanish in light of the trivial identity 




We see that the virial ground-state pressure (H^ agrees 
with that of the pressure obtained via the energy per 
particle [cf. Eq. (|4T]) ] provided that F{k) exists. Since 
the latter is a stronger condition than the existence of 
v{k), it is possible to devise a stealthy function v{k) for 
which F(k) does not exist and hence a virial pressure 
that either diverges or is nonconvergent. For example, 
this problem occurs for the power-law potential ()15|) with 
m = 0 (step function) for any dimension d. By contrast, 
the virial pressure is always well defined for the overlap 
potential d26l) in any dimension . This example serves 
to illustrate the mathematical subtleties that can arise 
because of the long-ranged nature of stealthy potentials 
in direct space. 


V. ENSEMBLE THEORY EOR STEALTHY 
DISORDERED GROUND STATES: EXACT 
INTEGRAL CONDITIONS ON THE PAIR 
STATISTICS 


Here, we derive some exact integral conditions that 
must be obeyed by both the pair correlation function 
52(r) and the structure factor S{k) for stealthy ground 
states that apply to general ensembles. These analytical 
relations can be profitably employed to test correspond¬ 
ing computer simulation results. 


A. General properties 

In any stealthy ground state, the structure factor at¬ 
tains its minimum value S{k) = 0 for 0 < k < K, and 
hence has the form 

Sik) = e{k-K)[l + Q{k)l (51) 

where 0(x) is the Heaviside step function defined by Eq. 
(fTTl) and Q{k) = S{k) — 1 is a function that obeys the 
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inequality Q{k) > —1. Therefore, from Eq. ((61), we have 
that the Fourier transform of the total correlation func¬ 
tion h{r) has the form 

ph{k)=f{k)+P{k), (52) 

where 

f{k) = -e{K-k) (53) 

and 

P{k) = Q{k- K)Q{k). (54) 

It is noteworthy that the function f{k) is identical to the 
Mayer-/ function for an equilibrium hard-sphere system 
in direct space. 

Taking the inverse Fourier transform of Eq. (1521) yields 
the direct-space total correlation function, given by 


ph{r) = /(r) -k P(r) 


(55) 


where 




d/2 


Jd/2{Kr) 


(56) 


and 


P{r) = 


> 


(27r)= Jk 
/ ^ \ <i/2 

\2Trr ) 


Jd/ 2 {Kr) - p, 


(57) 


general form (ED for h{k) into the definition of the total 
correlation function h{r) as obtained from Eq. (fl^ . and 
expanding h(r) in a Taylor series around r = 0 through 
second order in the radial distance r, yields 

h{r) = h{r = 0) + ^ ^ ’ (60) 

where 

poo 

h{r = 0) = -2dx + 2(fx k‘^-^Q{k)dk (61) 

Jk 

and the corresponding curvature is 

Therefore, from Eq. (EH, we see that the pair correla¬ 
tion function at the origin is given by 

nOO 

g 2 {r = 0) = 1 — 2dx + ^d'^X / k‘^~^Q{k)dk (63) 

Jk 

Since g2(?') must be non-negative for all r, we have the 
following integral condition on Q{k): 


2d\ 



k‘^-^Q{k)dk > 2dx-l. 


Hence, this integral must be positive for 



(64) 


(65) 


where the lower bound on P{r) indicated in Eq. (1571) 
follows from the fact that h{r) > —1 for all r for any 
point pattern. It trivially follows that since ph(k = 0) = 
p /jjd h{r)dr = —1, the volume integral of P{r) must be 
zero, i.e.. 


j P{r)dr = 0. 
Js.<‘ 


(58) 


Less trivially, because the product v{k)P{k) is zero for all 
fc, by Parseval’s theorem, we have the integral condition 



v{r)P{r)dr = 0. 


(59) 


Thus, the functions v{r) and P(r) are orthogonal to one 
another. The exact integral conditions (155)) and (1591) can 
be used to test the accuracy of numerical methods that 
yield estimates of the pair correlation function. 


We also conclude from Eq. El that for h{r) or g 2 {r) 
to have positive curvature at the origin, Q{k) must obey 
the additional integral condition; 

poo 

{d + 2) k‘‘+^Q{k)dk < 1 . ( 66 ) 

Jk 

Finally, we note that when g 2 {r = 0) = 0, the results 
above yield the equality 

k‘^-^Qik)dk = ^^^ (67) 

and, because the curvature must be positive in this in¬ 
stance, the inequality (1551) must generally be obeyed. 

The inequality (1541) . conditional inequality (1551) . and 
conditional equality (1671) provide integral conditions to 
test the accuracy of numerical methods that yield esti¬ 
mates of the structure factor. 


B. Behavior of the pair correlation function near 
the origin 

It is instructive to determine the behavior of the pair 
correlation function 32 (j') for small r. Substitution of the 


VI. EXISTENCE OF STEALTHY DISORDERED 
DEGENERATE GROUND STATES 

It is noteworthy that any periodic crystal with a fi¬ 
nite basis is a stealthy ground state for all positive y up 
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to its corresponding maximum value Xmax (or minimum 
value of the number density Pmin) determined by its first 
positive Bragg peak hsragg [minimal positive wave vec¬ 
tor for which S'(k) is positive]. Tables TIV list the pair 
Xmax,Pmin for some common periodic patterns in one, 
two, three, and four dimensions, respectively, all of which 
are part of the ground-state manifold; see Appendix A 
for mathematical definitions. (The crystals denoted by 
Dia^ and Kag^^ are d-dimensional generalizations of the 
diamond and kagome crystals, respectively, for d > 2 
[^.1 While the mere existence of such periodic ground 
states does not provide any clues about their occurrence 
probability in some ensemble, we will use these results 
here to show how disordered degenerate ground states 


arise as part of the ground-state manifold for sufficiently 
small X- 

TABLE I. Maximum values of x and corresponding minimum 
values of p for certain periodic stealthy ground states in R 
with A = 1. The configuration with the largest possible value 
of Xmax (smallest possible value of pmin) corresponds to the 
integer lattice. 


Structure 

\max 

Pmin 

Integer lattice (Z) 

1 

^ =0.15915... 

Periodic with n-particle basis 

1 

n 

^ = (0.15915...)n 


TABLE 11. Maximum values of x and corresponding minimum values of p for certain periodic stealthy ground states in R^ 
with K = 1. The configuration with the largest possible value of Xmax (smallest possible value of pmin) corresponds to the 
triangular lattice. 


Structure 

\max 

Pmin 

Kagome crystal (Kagj) 
Honeycomb crystal (Dia 2 ) 
Square lattice (Z^ = Z^) 
Triangular lattice (A 2 = A 2 ) 

3^= 0-3022... 
21^= 0-4534... 
f = 0.7853 ... 
7^= 0-9068... 

= 0.06581... 
^ = 0.04387 ... 

47r-^ 

A = 0.02533 ... 

47r-^ 

^ = 0.02193... 


TABLE III. Maximum values of x and corresponding minimum values of p for certain periodic stealthy ground states in R® 
with K = 1. Here MCC refers to the mean-centered cuboidal lattice, which is a Bravais lattice intermediate between the BCC 
and FCC lattices, and has an equivalent dual lattice [b^. The configuration with the largest possible value of Xmax (smallest 
possible value of pmin) corresponds to the BCC lattice. 


Structure 

Pyrochlore crystal (Kagg) 
Diamond crystal (Dias) 
Simple hexagonal lattice 
SC lattice (Z 3 = Z^) 
HCP crystal 
FCC lattice (D 3 = A 3 ) 
MCC lattice 
BCC lattice (D 3 = D 3 ) 


4\/T2 

TT 

2v/T2 

TT _ 


= 0.2267. 


= 0.4534. 


9 

2n 


0.6045.. 


= 0.6981... 


Sy/e-n- 


81 


0.7600 . 


VT2 


= 0.9068. 


0.9258... 


2V27r 


0.9873 . 


Pmin 


SvAir 

_1_ 

3V3n 


4^3 71 


= 0.01241... 
= 0.00620 ... 

■ = 0.00465 ... 


= 0.00403 ... 


3 ^= 0 - 00370 ... 
T = 0.00310 ... 


6^/37, 




0.00303... 

= 0.00285 ... 


TABLE IV. Maximum values of x and corresponding mini¬ 
mum values of p for certain periodic stealthy ground states in 
R"^ with K = 1. The configuration with the largest possible 
value of Xmax (smallest possible value of pmin) corresponds 
to the four-dimensional checkerboard lattice D 4 , = D%. 


Structure 

^max 

Pmin 

Kag 4 crystal 

= 0.2467... 

Ar = 0.001640 ... 

Dia 4 crystal 

^ =0.6168... 

-A- = 0.0006416 ... 

iDTr*^ 

Z"^ lattice 

=0.6168... 

jA = 0.0006416 ... 

D 4 lattice 

^ = 1.2337... 

3 A = 0.0003208 ... 


At fixed d, the smallest value of Pmin listed in Tables 
I-IV, which we call p^m^ corresponds to the dual of the 
densest Bravais lattice in direct space, and represents 
the critical density value below which a stealthy ground 
state does not exist for all k < \^Bragg\- The fact that 
Pmin corresponds to the body-centered-cubic (BCC) lat¬ 
tice for d = 3 was initially shown analytically in Ref. [l^ 
and subsequently numerically in Ref. 0- We note that 
the values of Pmin for the simple hexagonal lattice and 
hexagonal close-packed crystal for d = 3 reported in Ref. 
0 are incorrect because those calculations were based 
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on the erroneous assumption that the structure factors at 
the corresponding shortest reciprocal lattice vectors have 
nonvanishing values. 

Observe that in the case d = 1, there is no non-Bravais 
lattice (periodic structure with a basis n > 2) for which 
Xmax is greater than 1 / 2 , implying that the ground-state 
manifold is nondegenerate (uniquely the integer lattice) 
for 1/2 < y < 1. This case is to be contrasted with the 
cases d > 2 where the ground-state manifold must be 
degenerate for 1/2 < y < Xmax nondegenerate 
only at the point y = Xmax^ implied by Tables II-IV. 
Here, Xmax largest possible value of Xmax in some 

fixed dimension. 

Lemma .—At fixed K, a configuration comprised of the 
union (superposition) of m different stealthy ground- 
state configurations in with yi, y 2 , ■ ■ •, Xmi respec¬ 
tively, is itself stealthy with a y value given by 


retically. In particular, our objective is to derive ana¬ 
lytical formulas for the pair statistics of stealthy disor¬ 
dered ground states for sufficiently small y in the canon¬ 
ical ensemble as temperature T tends to zero; i.e., the 
probability of observing a configuration is proportional 
to exp[—d>(r'^)/(fcBr)] in the limit T ^ 0. We show 
here that under such circumstances, the pair statistics 
in the thermodynamic limit can be derived under the 
ansatz that stealthy ground states behave remarkably 
like pseudo-equilibrium hard-sphere systems in Fourier 
space. This ansatz enables us to exploit well-known ac¬ 
curate expressions for the pair statistics in direct space. 
As will be shown, agreement with computer simulations 
is excellent for sufficiently small y. 


A. “Pseudo” hard-sphere ansatz 


-1 

( 68 ) 

which is the harmonic mean of the Xi divided by m. 

Proof .—Formula (l68l) is a direct consequence of the 
fact that y is inversely proportional to the number den¬ 
sity p = Pi of union of the configurations in 
where pi is the number density associated with the ith 
configuration, which is inversely proportional to Xi- 
This Lemma, together with the fact that any periodic 
crystal with a finite basis is a stealthy ground state can 
be used to demonstrate rigorously how complex aperi¬ 
odic patterns can be ground states, entropically favored 
or not. A sketch of such a proof would involve the con¬ 
sideration of the union of m different periodic structures 
in with densities pi, p 2 , ■ • ■ j Pm, respectively, each of 
which are randomly translated and oriented with respect 
to some coordinate system such that m is very large but 
bounded and pi ^ pj for all i and j. It is clear that the 
resulting configuration will be a highly complex aperiodic 
structure in R'^ that tends toward a disordered stealthy 
pattern with a value of y that is very small but positive 
according to relation (1551) . 


X = 


Ha. ' 


VII. PAIR STATISTICS IN THE CANONICAL 
ENSEMBLE: “PSEUDO” HARD SPHERES IN 
FOURIER SPACE 


The task of formulating an ensemble theory that yields 
analytical expressions for the pair statistics of stealthy 
degenerate ground states is highly nontrivial because the 
dimensionality of the conhguration space depends on the 
density (or y) and there is a multitude of ways of sam¬ 
pling the ground-state manifold, each with its own prob¬ 
ability measure for finding a particular ground-state con¬ 
figuration. Therefore, it is desirable to specialize to equi¬ 
librium ensembles with Gibbs measures because the char¬ 
acterization of the ground states (as well as the corre¬ 
sponding excited states) would be most tractable theo¬ 


We have already noted that the step-function contribu¬ 
tion to ph{k) for stealthy ground states, denoted by f{k) 
in relation (l52l) . is identical to the Mayer-/ function for 
an equilibrium hard-sphere system in direct space. This 
implies that the corresponding contribution to S{k) is a 
simple hard-core step function 0(/c — K), which can be 
viewed as an equilibrium hard-sphere system in Fourier 
space with “spheres” of diameter K in the limit that y 
tends to zero. Why is this the case? Because such a 
step function is exactly the same as the pair correlation 
g 2 {r = k) of an equilibrium hard-sphere system in direct 
space in the limit that p tends to zero. That the structure 
factor must have the behavior S{k) —J- 0(fc — K) in the 
limit y —>■ 0 is perfectly reasonable, since a perturbation 
about the ideal-gas limit [where S{k) = I for all k] in 
which an infinitesimal fraction of the degrees of freedom 
are constrained should only introduce an infinitesimal 
change in S{k) of zero inside the exclusion zone (con¬ 
strained region). We call this the weakly constrained 
limit, where a step function S{k) is expected on maxi¬ 
mum entropy grounds; it corresponds to the most disor¬ 
dered (decorrelated) form of S{k) subject to the impene¬ 
trability condition in Fourier space. We refer to this phe¬ 
nomenon as equilibrated pseudo hard spheres in Fourier 
space because there are actually no points in that space 
that have a hard-core repulsion like true hard spheres do 
in direct space. 

On the same maximum entropy grounds, we expect 
that a perturbation expansion about the weakly con¬ 
strained limit y = 0 will lead to a perturbation expansion 
in y for S{k) that can be mapped to the low-density ex¬ 
pansion of 52 (c) for equilibrium hard spheres. More gen¬ 
erally, we make the ansatz that, in the canonical ensemble 
as T —s- 0, this hard-sphere analogy continues to hold as y 
is increased from zero to positive values, provided that y 
is small enough, implying that the collective coordinate 
variables n{k) (defined in the Introduction) are weakly 
correlated. Though the pseudo-hard-sphere picture must 
break down in some intermediate range of y, for d = 1 
and d = 2, this hard-sphere mapping is again exact when 
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X = Xmax, which corresponds to the maximal value of 
the packing fraction rj in these dimensions (see Tables 
I and II). This exact correspondence with the maximal 
value of r] when x = Xmax does not hold for d = 3 or 
d = 4, however. Thus, one should only expect that x 
and r] are proportional to one another, even at small x 
values. 

Under the pseudo-hard-sphere ansatz, the direct-space 
pair correlation function g 2 ^{r;r]) of a disordered hard- 
sphere system at a packing fraction rj for sufficiently small 
r] can be mapped into the structure factor S{k;x) for a 
disordered stealthy ground state derived from the canon¬ 
ical ensemble at fixed x for sufficiently small x as follows: 

S{k-,x) = 9?^{r = k]ri). (69) 

As alluded to above, the parameter x can be viewed as 
an effective packing fraction for pseudo hard spheres of 
diameter K in reciprocal space that is proportional to rj, 
i.e., 


V = Kd)x, (70) 

where 6(d) is a d-dependent parameter that is to be de¬ 
termined. Let hHs{f) be the total correlation function 
of a disordered equilibrium hard-sphere system in direct 
space and let us define for stealthy ground states 

H{k) = S{k) - 1 = ph{k). (71) 

The ansatz is also defined by the alternative mapping 

H(k) = hHs{r = k). (72) 

This mapping then enables us to exploit the well-known 
statistical-mechanical theory of equilibrium hard-sphere 
systems. In particular, we can employ a generalized 
Ornstein-Zernike convolution relation that defines the 
appropriate direct correlation function C{k), namely, 

H{k) = C{k) + r]H{k)®C{k), (73) 


where the symbol ® denotes the convolution operation 
in K.'^. Therefore, in direct space, H{r) is given by the 
relation of the following form: 


H{r) 


C{r) 

1 — (27r)^ rjC{r) 


For example, for d = 1, 




(74) 

(75) 


Inverting this function yields 


is highly accurate for low to intermediate densities along 
the liquid branch, or when mapped to the stealthy prob¬ 
lem, for low to intermediate values of x- 

It is noteworthy that the exact low-density expansion 
of hnsi’’’), foi' practical purposes, is sufficient to produce 
accurate estimates of H{k) = ph[k) and its counterpart 
ph{r) for low to intermediate values of x or rj. In particu¬ 
lar, using the mapping ((72ll . we obtain, for any dimension 
d, the following low-y expansion of ph(k) : 

ph{k) = -Q{K -k)[l + b{d) a{k; K)x + 0(x')] , 

(77) 

where 6(d) is the proportionality constant in Eq. (1701) and 
a{k; K) is the scaled intersection volume of two identical 
d-dimensional spheres of diameter K whose centers are 
separated by a distance fc [cf. Eq. (1001) ] This formula 
indicates that S{k) develops a peak value at fc = AT (over 
and above the value of unity due to the step function in 
the limit X —>■ 0) and then monotonically decreases until 
k = 2K, where it achieves its long-range value of unity 
for all k > 2K, which we will see is verified by computer 
simulations. Fourier inversion of Eq. dzzl), division by p, 
and use of (1351) yields a corresponding low-y expansion of 
the total correlation function h{r) through second order 
in X a-nd hence has an error term of order 

To get an idea of the large-r asymptotic behavior of 
the pair correlations, consider the limit x 0 for any d. 
In this limit, the total correlation function h{r) for any 
r obtained from Eq. dZZl) is given by 

PHr) = - Jd/ 2 {Kr) (x 0), (78) 

which for large r is given asymptotically by 

- ^(d+i )/2 -{d+ l)7r/4) (r^-Poo). 

(79) 

Thus, the longed-ranged oscillations of h[r) are con¬ 
trolled by the power law —Equation (1751) indi¬ 
cates that in the limits x 0 and p oo, h{r) —>■ 0, and 
therefore, the pair correlation function tends to the ideal 
gas even though the structure factor [Eq. (15^ ] cannot 
tend to the ideal-gas form because of its stealthy prop¬ 
erty. This result is in contrast to the situation considered 
in Fig. [1] where we take the p ^ oo limit by fixing N 
and letting vp 0. In that case, both 52 (j') and S{k) 
tend to the associated ideal-gas forms, i.e., g 2 {r) = 1 for 
all r and S{k) = 1 for all k. 

B. Comparison of theoretical predictions to 
simulations 


-rsin(r) -f (r[sin(r) -f cos(r)] - 1)?7 

- <™> 

For d = 2 and d = 3, one can use the Percus-Yevick clo¬ 
sure of the Ornstein-Zernike integral equation (s^ j. which 


In order to test our theoretical results for the pair 
statistics of stealthy ground states in the canonical en¬ 
semble, we have carried out computer simulations to gen¬ 
erate and sample such configurations, the details of which 
are described in Appendix B. In all cases, we take K = 1, 
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which sets the length scale. Our simulation results reveal 
that the functional trends for S{k) and g 2 i'r) predicted 
by the ansatz of pseudo hard spheres in Fourier space 
with an effective packing fraction x ^re remarkably ac¬ 
curate for a moderate range of x about x = 0- Because 
it is theoretically highly challenging to ascertain the pro¬ 
portionality constant b{d) in Eq. (1701) that arises in Eq. 
m, we must rely on the simulations to guide us in its 


determination. First, we observe that for d = 1, the 
mapping between x and r] is one to one, i.e., 6(1) = 1. 
Second, the simulation data suggest that, to an excellent 
approximation, b{d) for d > 2 is given by assuming that 
the peak value of S{k) or ph{k), achieved at k = K for 
sufficiently small x, is invariant with respect to this peak 
value as in the one-dimensional case, and consequently 
bid) = [a{K-K)2‘^]-\ 



C/D 



FIG. 6. Comparison of theoretical and simulation results for the structure factor S{k) for x = 0.05, 0.1, and 0.143 for d = 3. 
Here, K = 1. 




FIG. 7. Comparison of theoretical and simulation results for the pair correlation function g 2 (r) for x = 0.05, 0.1, and 0.143 
across the first three space dimensions. Here, K — 1. 


Figure [6] shows that the structure factor S{k), as ob¬ 
tained from Eqs. (El) and (ED, is in excellent agreement 
with the corresponding simulated quantities for x = 0.05, 
0.1, and 0.143 for d = 3. In Fig. [71 we compare our the¬ 
oretical results for the pair correlation function g 2 if), as 
obtained by Fourier inversion of Eq. (ED, to correspond¬ 
ing simulation results across the first three space dimen¬ 
sions. Again, we see excellent agreement between theory 
and simulations, which validates the pseudo-hard-sphere 
Fourier-space ansatz. Figure [8] depicts our theoretical 
predictions for g 2 i'r) for x = 0-15 across the first four 
space dimensions. It is seen that increasing dimensional¬ 
ity increases short-range correlations. 



FIG. 8. Theoretical predictions for the pair correlation func¬ 
tion 52 (r) for X = 0.15 across the first four space dimensions. 
Here K = 1. 


C. Translational order or disorder metric 

We have seen that both short- and long-scale corre¬ 
lations increase as x increases. A useful scalar positive 
order metric that captures the degree to which transla- 
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tional order increases with x is given by 

I/'<*“> 

where we have used Parseval’s theorem and D is some 
characteristic length scale [d^. Note that for an ideal 
gas (spatially uncorrelated Poisson point process), r = 0 
because h(r) = 0 for all r. Thus, a deviation of r from 
zero measures translational order with respect to the fully 
uncorrelated case. Because r diverges for any perfect 
crystal, it is a quantity that is better suited to distinguish 
the degree of pair correlations in amorphous systems. 

In the case of stealthy ground-state configurations, r 
is given explicitly by the relation 


where H{k) is given by Eq. (|7T|) . Substitution of the 
leading-order term in the x expansion Eq. dzii) into (HID 
yields 


T = 


4^2 (27r)'^ 


x" + 0{x^), 


(82) 


where we have taken D = K~^. Thus, for stealthy 
ground states, the order metric r grows quadratically 
with X for small x- Since the error is of order we 
expect that this quadratic form will be a very good ap¬ 
proximation of r up to moderately large values of x- Ih" 
deed, this is confirmed by our simulations up to x = 0.25. 
Note that because stealthy disordered ground states (for 
sufficiently small x) £^re pseudo-equilibrium hard-sphere 
systems in Eourier space, the form of the order metric r 
[Eq. (HU] ensures that it will behave similarly to r for 
equilibrium hard spheres in direct space for low densities. 


where vi{R) is the d-dimensional volume of a spherical 
window [cf. Eq. (1361) ]. h{r) is the total correlation func¬ 
tion [cf. Eq. ([U], a{r\R) is the scaled intersection vol¬ 
ume of two spherical windows of radius i?, as given by 
Eq. (l20l) . and d(fc; R) is the Eourier transform of a(r; R), 
which is explicitly given by 


d(fc;i?) =2‘^7r'^/2r(l-hd/2)-t^^^^^^^^. (84) 

rv 


We have already noted that the stealthy ground states 
considered in the present paper are hyperuniform, i.e., 
S'(fc) —)> 0 as fc —>■ 0 (see Sec. Ill Al) . This means that 
such systems obey the sum rule p h{r)dr = —1 and, 
because of the rapid manner in which S{k) vanishes in 
the limit fc —>■ 0 , the number variance has the following 
large-i? asymptotic behavior : 

a^{R) = A{R)R‘^-^ + OiR‘^-^), (85) 

where A{R) is a bounded function that oscillates around 
an average value 


A = lim — f A{R)dR. ( 86 ) 

L->oo L Jq 

The scaling (1851) occurs for a broader class of hyperuni¬ 
form systems, as specified by relation (fTOl) . The param¬ 
eter A is an order metric that quantifies the extent to 
which large-scale density fluctuations are suppressed in 
such hyperuniform systems [s^. To compare different 
hyperuniform systems, Torquato and Stillinger used the 
following rescaled order metric: 


B = 


A 


(87) 


VIII. LOCAL NUMBER VARIANCE FOR 
STEALTHY HYPERUNIFORM DISORDERED 
GROUND STATES 


Here we investigate theoretically the local number vari¬ 
ance for stealthy disordered ground states as a function of 
X and then use these results to extract an order metric 
[3^ that describes the extent to which large-scale den¬ 
sity fluctuations are suppressed as x increases in these 
hyperuniform systems (see Sec. Ill Al) . The local number 
variance cr^ (i?) associated with a general statistically ho¬ 
mogeneous and isotropic point process in at number 
density p for a spherical window of radius R is determined 
entirely by pair correlations [3^ : 


a'^{R) = pvi{R) 1 + p h{r)a{r; R)dr 


f S{k)a{k;R)d'k , (83) 


which is independent of the density, where (f) = pvi{l/2). 
Among all hyperuniform point patterns having the scal¬ 
ing (1551) . B is minimized (greatest suppression of large- 
scale density fluctuations) for the integer, triangular, 
BCC, and D 4 lattices for d = 1, 2, 3, and 4, respectively 

[nisi. 

Using the analytical results for g 2 {r) or S{k) described 
in the previous section, we have computed relation ((83l) 
for <t^{R) versus R for selected values of x across the 
first three dimensions and compared them to our corre¬ 
sponding simulation results; see Fig. |9l The analytical 
and numerical results are in excellent agreement with one 
another. Tabled lists the order metric B for various val¬ 
ues of X across the first four dimensions for disordered 
stealthy ground states, as obtained from the analytical 
estimates of (j‘^{R) and Eq. (1571) . These results are also 
compared to the corresponding optimal values. As ex¬ 
pected, B decreases as x increases for fixed d. 
















16 




FIG. 9. Comparison of analytical and numerical results for 
the number variance for x = 0.05, 0.1, and 0.2 across the hrst 
three space dimensions. Here, K — 1. 


TABLE V. The order metric B for various values of x across 
the first four dimensions for disordered stealthy ground states 
for selected values of y up to x = 0.25, as obtained from 
the analytical estimates of ct^{R) and Eq. (1871) . Included for 
comparison are the structures in each dimension that have the 
minimal values of B, all of which are Bravais lattices [s^, [s^. 


X 

d= 1 

d= 2 

CO 

II 

0.05 

2.071 

1.452 

2.164 

4.560 

0.1 

1.051 

1.040 

1.738 

3.875 

0.143 

0.745 

0.880 

1.558 

3.576 

0.2 

0.54 

0.755 

1.411 

3.327 

0.25 

0.439 

0.683 

1.325 

3.179 

Integer lattice 
Triangular lattice 

BCC lattice 

U 4 lattice 

0.167 

0.508 

1.245 

2.798 


IX. NEAREST-NEIGHBOR FUNCTIONS 

Here, we obtain theoretical predictions for the nearest- 
neighbor functions of stealthy disordered ground states. 


Nearest-neighbor functions describe the probability of 
finding the nearest point of a point process in at some 
given distance from a reference point in space. Such sta¬ 
tistical quantities are called “void” or “particle” nearest- 
neighbor functions if the reference point is an arbitrary 
point of space or an actual point of the point process, re¬ 
spectively Our focus here is on the particle nearest- 
neighbor functions. 

The particle nearest-neighbor probability density func¬ 
tion Hp{r) is defined such that Hp{r)dr gives the prob¬ 
ability that the nearest point to the arbitrarily chosen 
point lies at a distance between r and r dr from this 
chosen point of the point process. The probability that 
a sphere of radius r centered at a point does not have 
other points, called the exclusion probability Ep{r), is 
the associated complementary cumulative distribution 
function and so Ep(r) = 1 — L Hp(x) dx and and hence 
Hp{r) = -dEpjdr. 

The nearest-neighbor functions can be expressed as an 
infinite series whose terms are integrals over n-body cor¬ 
relation functions defined in Sec. inilnl,[2l- In general, 
an exact evaluation of this infinite series is not possible 
because the are not known accurately for n > 3, ex¬ 
cept for simple cases, such as the Poisson point process. 
Theoretically, one must either devise approximations or 
rigorous bounds to estimate nearest-neighbor quantities 
for general models 0, [t^ . 

Torquato has given rigorous upper and lower bounds 
on the so-called canonical n-point correlation function 
Hn for point processes in R'* [tI, [tJ. Since nearest- 
neighbor functions are just special cases of Eln, then we 
also have strict bounds on them for such models 0,123. 
Here, we employ upper and lower bounds on Ep{r), 
which relies on knowledge of the pair correlation func¬ 
tion: 


Ep{r)>l-Z{r), (88) 

Ep{r) < exp[-Z(r)], (89) 


where Z{r) is the cumulative coordination number [cf. 
dH)]. The upper bound (IMl) was presented in Ref. 

These bounds are evaluated for stealthy ground-state 
configurations using the analytical expression for the pair 
correlation function given in Sec. IVHI Figures [TUI and [TT] 
compare these bounds to our numerical results for both 
X = 0.05 and y = 0.1 for d = 2 and d = 3, respectively. 
We see that the bounds on Ep{r) provide the correct 
qualitative trends as a function of r; the upper bound 
being the sharper of the two bounds for these cases. 
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FIG. 10. Comparison of the lower and upper bounds (I88II and (1891) to our numerical results for the exclusion probability 
function Ep{r) for x ~ 0-05 and x = 0-1 for d = 2. Here, K = 1. 




FIG. 11. Comparison of the lower and upper bounds (I88II and (1891) to our numerical results for the exclusion probability 
function Ep{r) for x = 0.05 and x = 0.1 for d = 3. Here, K = 1. 


The mean nearest-neighbor distance A is defined as the 
first moment of Hp{r) or, equivalently, zeroth moment 
of Ep{r), i.e.. 


A = / rHp{r)dr = / Ep{r)dr. (90) 
Jo Jo 


For an ideal gas (Poisson point process) at number den¬ 
sity p, the mean nearest neighbor can be explicitly given 
in any dimension [^ : 



FIG. 12. Upper bounds on the mean nearest-neighbor dis¬ 
tance A, scaled by Xueai [cf. Eq. (|9T1) ]. as a function of x for 
the first four space dimensions. 


^Ideal 


r(i + i/d) 

2[pui(l/2)]i/d- 


(91) 


Using the upper bound (l8^ and relation (1^ . we plot 
in Fig. [12] upper bounds on the mean nearest-neighbor 
distance A, scaled by the corresponding ideal-gas quan¬ 
tity obtained from Eq. (ED, as a function of x for the hrst 
four space dimensions. For fixed Xi tire upper bounds on 
^/Kdeai decrease as the space dimension increases, as ex¬ 
pected, and tends to unity in the large-d limit, consistent 
with the so-called “decorrelation” principle . 


X. EXCITED STATES: STRUCTURE FACTOR 
AND THERMAL EXPANSION COEFFICIENT 


Here, we derive accurate analytical formulas for the 
structure factor and thermal expansion coefficient for the 
excited states associated with stealthy ground states at 
sufficiently small temperatures. We see from the com¬ 
pressibility relation ([3D that if the isothermal compress¬ 
ibility kt is bounded, then S'(O) must be zero for any 
ground state, stealthy or not. Recall that for stealthy 
ground states, up is bounded according to Eq. (1421) . 
Now consider excited states inhnitesimally close to the 
stealthy ground states, i.e., when temperature T is pos¬ 
itive and infinitesimally small. Under the highly plausi- 
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ble assumption that the structure of such excited states 
will be infinitesimally near the ground-state configura¬ 
tions for sufficiently small x and T, then to an excellent 
approximation, the pressure is given by 

P-pT+^, (92) 

where the first term is the ideal-gas contribution and the 
second term is the configurational contribution, which, 
under the stated conditions, is effectively the same as 
the ground-state expression dUD, where we have set 
ks = vq = K = 1. Thus, relation (IMl) yields the isother¬ 
mal compressibility kt = [p{p + T)]~"^, which, when sub¬ 
stituted into Eq. dS]) for large p (small y) and small T, 
yields that iS'(O) varies linearly with T for such excited 
states: 

S{0)^C{d)xT, (93) 

in units ks = vq = K = 1, where C{d) = 2d (27r)'^/wi(l) 
is a d-dependent constant. 

Figure [Ql shows that the prediction of relation (IMl) is 
in excellent agreement with our MD simulation results 
(Appendix B) in the case d = 2. It is expected that this 
positive value of S'(O) will be the uniform value of S{k) 
for 0 < fc < AT for the special case of the step-function 
power-law potential v{K) [the case m = 0 in Eq. dTSl) ] for 
small X- This behavior of S{k) has indeed been verified 
by our simulation results in various dimensions. For other 
stealthy-potential function choices, S{k) will no longer be 
a constant for 0 < k < K. 



FIG. 13. Comparison of the theoretically predicted structure 
factor at the origin S'(O) versus absolute temperature T, as 
obtained from Eq. (I93II . to our corresponding simulation re¬ 
sults for a two-dimensional stealthy system at y = 0.1. Here 
we take ks ~ vq = K = 1. 


An interesting conclusion to be drawn from this anal¬ 
ysis is that, for a system, S'(O) can be arbitrarily close 
to zero at positive temperatures, even if T is itself arbi¬ 
trarily small. This means that, for all practical purposes, 
such systems at positive T are effectively hyperuniform. 
Perfect hyperuniformity is not necessarily required in or¬ 
der to achieve novel physical properties in technological 
applications. 

Using the cyclic identity 

{dp/dT)p{dT/dp)p{dp/dp)T = —I and approxima¬ 
tion (1^^ . it immediately follows that the thermal 


expansion coefficient a = —p ^{dp/dT)p, for sufficiently 
small X and T, is given by 

a^C\d)x^T. ( 94 ) 

We see that the thermal expansion is positive under such 
conditions, which is to be contrasted with the anoma¬ 
lous negative thermal expansion behavior for sufficiently 
large y over a low temperature range demonstrated in 
our earlier numerical work [isl . 

XI. CONCLUSIONS AND DISCUSSION 

Stealthy hyperuniform disordered ground states in 
are infinitely degenerate and arise from a class of bounded 
long-ranged pair potentials with compactly supported 
Fourier transforms. Such exotic many-particle states of 
matter were previously studied only numerically. Be¬ 
cause the configurational dimensionality depends on the 
density (or y), a highly unusual situation, and there are 
an infinite number of distinct ways to sample the ground- 
state manifold, each with its own probability measure, 
it has been theoretically very challenging to devise pre¬ 
dictive ensemble theories. A new type of statistical- 
mechanical theory needed to be invented. This paper 
has initiated such a theoretical program. 

Specifically, we have derived general exact relations 
for the ground-state energy, pressure, and isothermal 
compressibility that apply to any ensemble as a func¬ 
tion of the number density p in any dimension d. We 
demonstrated how disordered degenerate ground states 
can arise as part of the ground-state manifold. We also 
obtained exact integral conditions that both the pair cor¬ 
relation function g2ir) and structure factor S{k) must 
satisfy in any ensemble. Then, we specialized our results 
to the canonical ensemble in the zero-temperature limit 
by exploiting an ansatz that stealthy states behave like 

H udo-equilibrium hard-sphere systems in Fourier space 
. The resulting theoretical predictions for g 2 ir) and 
S{k) were shown to be in excellent agreement with com¬ 
puter simulations across the first three space dimensions 
for sufficiently small y. These results were used to the¬ 
oretically obtain order metrics, local number variance, 
and nearest-neighbor functions across dimensions. We 
also derived accurate analytical formulas for the struc¬ 
ture factor and thermal expansion coefficient for the ex¬ 
cited states associated with stealthy ground states at suf¬ 
ficiently small temperatures. Our analyses provide new 
insights on our fundamental understanding of the nature 
and formation of low-temperature states of amorphous 
matter. Our work also offers challenges to experimental¬ 
ists to synthesize stealthy ground states at the molecular 
level, perhaps with polymers, as suggested in Sec. IIIIl 
There are many remaining open theoretical problems. 
While the pseudo-hard-sphere system picture for the 
canonical ensemble is almost surely exact in the limit 
y —>■ 0, a future challenge would be to provide rigorous 
justification for this picture for positive but small y. One 
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possible avenue that could be pursued is the formulation 
of an exact perturbation theory for the pair statistics 
about the weakly constrained limit (y 0). As noted 
in the Introduction, while the configuration space is fully 
connected for sufficiently small y, quantifying its topol¬ 
ogy as a function of y up to Xmax outstanding open 
problem. At some intermediate range of y, the topol¬ 
ogy of the ground-state manifold undergoes a sequence 
of one or more disconnection events, but this process is 
poorly understood and demands future study. In the 
limit y —>■ Xmax^ the disconnection becomes complete at 
the unique crystal ground state [t^ . 


Disordered 


Crystal 


X=0(P=oo) X=l/2 XLx(P=Pmin) 


FIG. 14. Phase diagram for the entropically favored stealthy 
ground states in the canonical ensemble as a function of y 
which applies to the first four space dimensions. 


but this remains to be confirmed. 

All of our previous and current simulations for the first 
three space dimensions strongly suggest that all of 

the energy minima attained were global ones for y < 0.5, 
but when y > 0.5, the topography of the energy land¬ 
scape suddenly exhibits local minima above the ground- 
state energies. The possible configurations that can arise 
as part of the ground-state manifold for y > 1/2, regard¬ 
less of their probability of occurrence, not only include 
periodic crystals for d > 2, as discussed in Sec. rvn but 
also aperiodic structures, reflecting the complex nature 
of the energy landscape. For example, for d = 2, the 
manifold includes aperiodic “wavy phases,” which have 
been shown to arise via numerical energy minimizations 
from random initial conditions with high probability in a 
range of y wh ere the triangular lattice is entropically fa¬ 
vored [IJ, [Ij . A deeper understanding of such aspects of 
the ground-state manifold would undoubtedly shed light 
on the topography of the energy landscape. 


The simple constraint and degrees-of-freedom count¬ 
ing arguments described earlier lead to definite predic¬ 
tions for the entropically favored stealthy ground states 
derived from the canonical ensemble in the limit T ^ 0. 
For y between 0 and 1/2, the ground states are disor¬ 
dered and possess a configurational dimension per par¬ 
ticle dc of d{l — 2y); see Eq. (l37)) . At y = 1/2, the 
configurational dimensionality per particle collapses to 
zero, and there is a concomitant phase transition to a 
crystal phase. The fact that the pressure is a continuous 
function of density [cf. Eq. (|4T|)] for all y up to Xmax 
(see Tables I-IV) implies that any phase transition could 
be continuous. While this eliminates a first-order phase 
transition in which the phase densities are unequal, it 
does not prohibit such a phase transition in which the 
two distinct phases possess the same number density. At 
y = 0.5, our numerical evidence indicates that the struc¬ 
turally distinct fluid and crystal phases have equal free 
energies. This does not necessarily imply that the two 
phases could coexist side by side within the system sepa¬ 
rated by an interface. This phase diagram is depicted in 
Fig. m which applies to the first four space dimensions. 
For d = 1, the only crystal phase allowed is the integer 
lattice (see Sec. ED, and hence there can be no phase 
coexistence. For d = 2, our simulations |14l - [l6l [T^ in¬ 
dicate that the crystal phase is the triangular lattice for 
1/2 < y < Xmax- However, for d = 3, it is possible that 
there may be more than one crystal phase. For exam¬ 
ple, while for sufficiently high y up to Xmax^ expect 
the stable crystal to be the BCC lattice, our simulations 
cannot eliminate the possibility that the FCC lattice is a 
stable phase for some y in the range 1/2 < y < 0.9068 ...; 
for y > 0.9068..., the FCC lattice cannot be a ground 
state (see Table III). Since four dimensions is more sim¬ 
ilar to two dimensions in that the lattice corresponding 
to Xmax is equivalent to its dual, we would expect that 
the Du lattice is the stable crystal for 1/2 < y < Xmax^ 
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Appendix A: COMMON d-DIMENSIONAL 
LATTICES 

Common d-dimensional lattices include the hypercu- 
bic checkerboard Dd, and root Ad lattices, defined, 
respectively, by 

= {{xi,... ,Xd) ■ Xi gZ} for d > I (Al) 

Dd = {{xi,..., Xd) G Z‘^ : xi + ■ ■ ■ + Xd even} for d > 3 

(A2) 

Ad = {(a;o, xi,... ,Xd) G Z'^+^ : xq + xi -I-h Xd = 0} 

for d > I, (A3) 

where Z is the set of integers (.. . — 3, —2, —1, 0,1, 2,3 ...); 
xi,... ,Xd denote the components of a lattice vector of 
either Z'^ or Dd; and xg, xi,..., Xd denote a lattice vector 
of Ad. The d-dimensional lattices Z/, and A^ are 
the corresponding dual lattices. Following Conway and 
Sloane [Tg , we say that two lattices are equivalent or 
similar if one becomes identical to the other possibly by a 
rotation, reflection, and change of scale, for which we use 
the symbol =. The Ad and Dd lattices can be regarded as 
d-dimensional generalizations of the face-centered-cubic 
(FCC) lattice defined by Ag = Dg-, however, for d > 4, 
they are no longer equivalent. In two dimensions, A 2 = 
A 2 defines the triangular lattice with a dual lattice that 
is equivalent. In three dimensions, A| = Dg defines the 
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body-centered-cubic (BCC) lattice. In four dimensions, 
the checkerboard lattice and its dual are equivalent, i.e., 
ZI 4 = D|. The hypercubic lattice = Z^ and its dual 
lattice are equivalent for all d. 


Appendix B: SIMULATION PROCEDURE TO 
GENERATE AND SAMPLE STEALTHY 
GROUND STATES IN THE CANONICAL 
ENSEMBLE 

To numerically sample the stealthy ground-state mani¬ 
fold in the disordered regime in the canonical ensemble in 
the T 0 limit for d = 1, 2, and 3, we performed molec¬ 
ular dynamics (MD) simulations at a very low dimen¬ 
sionless equilibration temperature Te = kBT/[voK^), 
periodically took configurational “snapshots,” and then 
used these configurations as input to the L-BFGS opti¬ 
mization algorithm to get the corresponding ground 
states. The dimensionless temperatures that we use are 
Te = 2 X 10-4 for d = 1, Te = 2 X 10"® for d = 2, and 
Te = 1 X 10“® for d = 3. The equilibration tempera¬ 
ture at a fixed dimension was chosen so that no changes 
in the pair correlation function are observed over some 
range of equilibration temperatures. The MD simula¬ 
tions were first performed in the microcanonical ensemble 
using the velocity Verlet algorithm [s^. The time steps 
were chosen so that the relative energy change every 3000 
time steps is less than 10“®. However, to enforce the de¬ 
sired temperature, we also performed MD simulations in 
the canonical ensemble using an Anderson thermostat 


[80| . We employed fifteen million time steps to equili¬ 
brate a system. After that, a snapshot was taken every 
3000 time steps for further energy minimization. Because 
stealthy potentials in direct space are long-ranged, the 
energy is most accurately calculated in Fourier space us¬ 
ing Eq. (1). The numerical errors in the achieved ground- 
state energies are extremely small, usually on the order 
of 10“^® (in units of vqK'^). The force on the jth parti¬ 
cle is calculated using the gradient of Eq. (1), yielding 
^3 = = :;^Ekki*(k)Im[h(k)exp(ik-rj)]. 

The number of particles in the simulation box, TV, is cal¬ 
culated from Eq. (2) for given d, x and M{K). We 
chose M{K) = 50 for d = 1, M{K) = 54 for d = 2, and 
M{K) = 39 for d = 3. The fundamental cell employed 
is the one corresponding to the crystal with the largest 
value of Xmax in each dimension (see Tables I-IV). 

Structural characteristics, such as the pair correlation 
function 52 (f), structure factor S{k), number variance 
(j‘^{R) and nearest-neighbor function Ep(r), by sampling 
each generated configuration for a fixed value of x and en¬ 
semble averaging over at least 20000 conhgurations. Two 
power-law potentials (13 were used: one with m = 0 and 
the other with m = 2 . As expected, both potentials pro¬ 
duced the same ensemble-averaged structural properties 
to within small numerical errors (as explained in Sec. 
IIVB|) . the agreement of which provides a good test on 
the validity of the simulation results. Our simulation re¬ 
sults for g 2 {r) and S{k) also satished the exact integral 
conditions presented in Secs. IV Al and IV Bl Additional 
simulation details will be described elsewhere 0 . 
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